Probing electronic excitations in molecular conduction 
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We identify experimental signatures in the current-voltage (I-V) characteristics of weakly con- 
tacted molecules directly arising from excitations in their many electron spectrum. The current is 
calculated using a multielectron master equation in the Fock space of an exact diagonalized model 
many-body Hamiltonian for a prototypical molecule. Using this approach, we explain several non- 
trivial features in frequently observed I-Vs in terms of a rich spectrum of excitations that may be 
hard to describe adequately with standard one-electron self-consistent field (SCF) theories. 
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Theoretical efforts to model molecular conduction have 
largely been based on SCF models for electron-electron 
interactions [l|. While they have been fairly successful 
in describing both shapes and magnitudes of various I- 
V characteristics 0, , notable exceptions include low- 
temperature measurements on unconjugated and weakly 
coupled molecules 0, 0, El EL as we ^ as snor t conju- 
gated molecules where there are clear disagreements 
between theory and experiment. Some disagreements 
could be attributed to uncertainties in geometry or par- 
asitic resistances; nevertheless the applicability of SCF 
approaches need to be scrutinized, especially in the weak 
coupling regime. Charging energies of short molecules 
(~ 3 eV for benzene) are often larger than their contact 
induced broadenings (< 0.2 eV for benzene dithiol on 
gold), while couplings between various molecular units 
(~ 2 eV for conjugated molecules, much less for non- 
conjugated species) can be tuned widely using synthetic 
chemistry. It is thus debatable whether a molecule acts 
as a quantum wire in the SCF regime, or as a quantum 
dot in the Coulomb Blockade (CB) regime. 

In this paper, we employ a multielectron master equa- 
tion 0, [H| in the Fock space of a prototypical molecular 
Hamiltonian to describe conduction through molecules 
with weak contact couplings or poor conjugation. A 
full many-body treatment of transport even with a small 
molecule, modeled simply as an array of quantum dots, 
yields many features with compelling similarities (Fig. 1) 
to relevant experiments jj, |5|, |6|. These features, how- 
ever, are quite difficult to obtain using a traditional non- 
equilibrium Green's function (NEGF) treatment of trans- 
port, being only perturbative in the interaction param- 
eter 11]. A spin restricted (RSCF) calculation (inset 
in Fig. lc) typically creates slow current onsets spread 
over several volts by Coulomb costs for adiabatic charg- 
ing. The high zero-bias conductances, in clear variance 
with experiments, could be removed by incorporating 
self-interaction correction for integer charge addition in 
the CB regime. However, crucial to experiments in this 
regime is the fact that the molecule can also execute 
transitions between various excited states of the neutral 
and singly charged species at no additional Coulomb cost. 
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FIG. 1: (color online) (a, b) Experimental 0, 13 and (c, d) 
theoretical I-Vs for a molecular ring weakly coupled with a 
backbone or with conducting electrodes. Many nontrivial, 
features such as low zero-bias conductance, sharp current on- 
set, and a subsequent quasilinear region spanning several volts 
with multiple closely spaced features (a-d) arise in our treat- 
ment of CB from excitations, but not even qualitatively in a 
spin restricted SCF (RSCF) treatment for the same param- 
eter set |13| (inset in c). For asymmetric contacts, there are 
additional features (b, d) including current step heights (as 
opposed to widths) that are asymmetric in bias, are modu- 
lated with a gate voltage 0, and reverse polarity for gate 
voltages on either side of the charge degeneracy point 
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making it possible to directly probe a rich spectrum of 
such transition levels within a small bias window. 

It seems difficult to capture this rich spectrum ade- 
quately within any SCF theory even with self interaction 
correction 0, Q, or effective one electron poten- 
tials especially under non-equilibrium conditions. A 
single spin-degenerate level (Fig. 2) illustrates the prob- 
lem. While the deficiencies of SCF (e.g. adiabatically 
smeared steps) are rectified with self-interaction correc- 
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FIG. 2: (color online) (a) Fock space, (b) equilibrium occu- 
pancy N — 11 and (c) nonequilibrium I-V of a spin degen- 
erate level e = leV with a single electron charging energy 
U = leV. A restricted SCF (RSCF) (pink) calculation shows 
fractional charge occupation and is inappropriate in the weak 
coupling limit. A Spin Unrestricted SCF (blue) describes inte- 
ger charge transfer and matches the many-body N — fi (black) 
plot; however, it yields equal current steps corresponding to 
sequential removal of two electrons, as opposed to a many- 
body calculation (shown in black) in which step heights are 
in the ratio of 2:1. Including correlations in SCF alters the 
current onsets and the plateau widths, but misses the essen- 
tial point that consecutive removal of spins need not carry 
equal current |^.H^. 



tions using a spin-unrestricted calculation for equilibrium 
properties such as N — fi (Fig. 2b), the same approach 
gives the wrong nonequilibrium properties such as current 
step heights. An unrestricted calculation yields equal 
step heights for each spin removal, while the exact result 
using rate equations predicts that the first step is two- 
thirds of the second (Fig. 2c), there being two ways to 
remove (add) the first spin for a filled (empty) level, but 
only one way to remove (add) the second one (Fig 2a). 

The discrepancy with SCF becomes more pronounced 
for multiple orbitals where a spin can be removed by 
one contact from the ground state and reinjected by the 
other into feasible excitations of the neutral and singly 
charged systems, creating additional features within the 
Coulomb Blockade plateaus. Such excitations, crucial to 
the experiments addressed here, arise non-perturbatively 
from our rate equations through exact diagonalization 
of the many-body Hamiltonian, going beyond orthodox 
Coulomb Blockade theory [HI due to size quantization 
and transitions among discrete many-body states. As the 
size of multielectron Fock space increases exponentially 
with the number of basis functions, we employ a mini- 
mal basis set in a reduced single-particle Hilbert space 
that captures the conjugation chemistry and electrostat- 
ics and yet allows exact diagonalization Quantita- 
tive justice to chemistry would possibly require looking 
at a reduced subset of excitations (partial configuration 



interaction |2(J) within a multiorbital description. Our 
aim is to solve the transport problem exactly for a simple 
system, rather than do an approximate SCF calculation 
of a more elaborate quantum chemical system . 

Approach. We start with a tight-binding model for 
benzene (one orbital per atom), with onsite, hopping, 
and Hubbard parameters that can be benchmarked with 
separate LDA calculations |22|. In contrast with single 
dot studies, long-ranged Coulomb terms (modeled with 
the Mataga-Nishimoto approximation 23]) and hopping 
are responsible for off-diagonal correlations in the charg- 
ing term of the molecular eigenspace. Exact diagonaliz- 
ing this Hamiltonian yields a large spectrum of closely 
spaced excitations in every charged configuration of the 
molecule. Using the equation of motion of the density 
matrix of the composite molecule and leads and assuming 
no molecule-lead correlations, one can derive [2 HH a 
simple master equation for the density-matrix of the sys- 
tem. Ignoring off-diagonal coherences, we are left with 
a master equation |23 in terms of the occupation prob- 
abilities P[* of each N electron many-body state \N,i) 
with total energy Ef. The master equation then in- 
volves transition rates R(N,i)^(N±i,j) between states dif- 
fering by a single electron, leading to a set of independent 
equations defined by the size of the Fock space Q 

dP N 

= ~ [ R (N,i)^(N±lJ)P l N ~ R(N±l,j)^(N,i)Pj N±1 ] 
3 

(1) 

along with the normalization equation J2 i N P^ = 1. 
For weakly coupled dispersionless contacts, parameter- 
ized using bare-electron tunneling rates j a , (a: left/right 
contact), we define rate constants 

= 7a\(N,i\ci\N-lJ)\ 2 
rj£ = ~f a \(N,i\c a \N + l,j)\ 2 , (2) 

c4, c a are the creation/annihilation operators for an elec- 
tron on the molecular end atom coupled with the corre- 
sponding electrode. The transition rates are given by 

R(N,i)^(N-l,j) = ^2 ^ijL t 1 - f( e ij T - Va)] 
a>=L,R 

R(N-l,j)^(N,i) = ^2 ^ijaf( e ij T ~^a)' (3) 
a=L,R 

for the removal levels (TV, i — > N — 1, j), and replacing 
(r — > a, / — > 1 — /) for the addition levels (TV, i — > TV + 
1, j). ji a are the contact electrochemical potentials, / 
is the corresponding Fermi function, with single particle 
removal and addition transport channels e^J r = — 

Ef~\ and e% a = Ef +1 - E? . Finally, the steady-state 
solution to Eq. is used to get the left terminal current 

j u e \t>l p n n L pN±i~ 

ij 

(4) 
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FIG. 3: (color online) Coulomb Blockade I-V features for a 
general molecular system, Threshold transition involving 1) 
only ground states. Here \Ef — eocfl — 10 meV: e^ r is accessed 
before e^ r (shown in the adjacent state transition diagram). 
I-V characteristics (shown black in (b)) then has a brief in- 
tervening plateau until an excitation is accessed. 2) Excited 
states: Conduction at threshold (\Ef — €oo\ = 30 meV) in- 
volves a transport channel involving excited states also (say 
efo 7 *) i.e., e^ r is accessed before e^ r . In this case (see text) 
current rise (shown blue in (c)) due to closely spaced excita- 
tions follows upon threshold. 



where states corresponding to a removal of electrons by 
the left electrode involve a negative sign. 

Results. We calculate current in a break-junction 
configuration with equal electrostatic coupling with the 
leads, setting /il,r — T e ^/2, and equal resistive 
couplings set by the ratio 7 = Jl/jr — 1? 7l = 0.6 meV. 
Coulomb Blockade with integer charge transfer manifests 
itself as a vanishing pre-threshold current followed by a 
step wise increase in current 0, IM El ^3 • The onset of 
conduction is established by the offset between the equi- 
librium Fermi energy Ep and the first accessible transi- 
tion energy (focussing on removal levels for concreteness, 
this corresponds to the transport channel marked Cqq in 
Figs. 3a and 3c). The onset can be varied by varying 
the gate voltage, thereby accounting for the variation in 
conductance gap with gate bias. 

The simplest impact of Coulomb Blockade on the I-Vs 
of short molecular wires is a clear suppression of zero- 
bias conductance, often seen experimentally H,[2(|. In- 
deed, a spin , unrestricted SCF with self-interaction cor- 
rections 0, can yield a Coulomb staircase with in- 
tervening plateaus through the Coulomb cost of adding 
or removing an electron to the molecular ground state. 
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FIG. 4: (a) State transition diagram showing various addition 
and removal pathways for asymmetric contacts (71, ^> 7^), 
including the possibility of populating higher excitations (b), 
say, via transport channel e^ r at threshold. For positive bias 
charge removal is the rate limiting process, while for negative 
bias addition dominates, accounting thus for the correspond- 
ing I-V asymmetry in (c). Progressive access of higher exci- 
tations also accounts for the observed gate modulation of the 
current steps, as shown in (d). 



However, integer charge transfer can also occur between 
various electronic excitations of the neutral and singly 
charged species at marginal correlation costs |2?|. The 
above fact leads to fine structure in the plateau regions 
0, 0, 0, 0, specifically, a quasilinear regime resulting 
from very closely spaced transport channels (e^J) via ex- 
citations. The crucial step is the access of the first excited 
state via channel e^ r , following which transport channels 
involving higher excitations are accessible in a very small 
bias window. The sequence of access of transport chan- 
nels upon bias, enumerated in the state transition dia- 
grams shown in Figs. 3a and 3b, determines the shape 
of the I-V. When the Fermi energy Ep lies closer to the 
threshold transport channel eQQ r (Fig. 3a), it takes an 
additional positive drain bias for the source to access the 
first excited state of the neutral system via the transition 
efo r , as shown in the state transition diagram in Fig. 3a. 
The I-V shows a sharp rise followed by a plateau (Fig. 
3c), as seen in various experiments [2^. However, when 
transport channels that involve low lying excitations such 
as e^Q r are closer to the Fermi energy E F than e^ r (Fig. 
3b), the excitations get populated by the left contact im- 
mediately when the right contact intersects the thresh- 
old channel e^, allowing for a simultaneous population 
of both the ground and first excited states via e^ r and 
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at threshold. Under these conditions the I-V shows a 
sharp onset followed immediately by a quasilinear regime 
(Fig. 3c) with no intervening plateaus, as observed fre- 
quently in I-Vs of molecules weakly coupled with a back- 
bone n,n@. 

The direct role of excitations in conduction becomes 
particularly striking under asymmetric coupling (7 = 



100, 7l = 0.6 meV) with contacts 0, In contrast 

to the SCF regime where unequal charging drags out the 
same level current over different voltage widths Q , in the 
CB regime the current step heights themselves are asym- 
metric at threshold (Fig. 4c). This asymmetry arises due 
to the difference in the number of pathways for removing 
or adding a spin, also taking into account the possible ex- 
citation channels between the neutral and singly charged 
species (Fig. 4 a, b). The number of such accessible 
excitations at threshold can be altered with an external 
gate bias, leading to a prominent gate modulation of the 
threshold current heights, over and above the modulation 
of the onset voltages and the conductance gap Q (Fig. 
4d). Furthermore, it is easy to show that the asymme- 
try will flip between gate voltages on either side of the 
charge degeneracy point, as is also observed experimen- 
tally 12]. While the qualitative features of our I-Vs are 
robust with respect to variation of our model parameters, 
details specific to experiments (e.g. onset voltages, po- 
larization asymmetries and temperature dependences 
@|) can vary and will be discussed in detail elsewhere [22^ . 
For instance, correlation alone cannot explain ultralow 
peak currents through a level since those depend only 
on contact couplings through the ratio 7l7,r/(7l + Jr)- 
This predicts a peak current ~ 3 fiA for a O.leV broad- 
ening as in chemisorbed benzene dithiol [2(J, still much 
larger than experiments [g , indicating that one needs fur- 



ther to postulate weak couplings due to non-ideal contact 
couplings or perhaps parasitic resistances from multiple 
molecules [2<j. Further complications could arise from 
strong electron-phonon interactions Q that smoothen 
out the first few conduction plateaus (Fig. Qi) due to 
low lying excitations over phonon energies significantly 
smaller than their Coulomb counterparts at tens of meV. 

In summary, we have used a rate equation in the Fock 
space of a molecular Hamiltonian to address significant 
experimental features like suppressed zero-bias conduc- 
tances, sharp steps that are often asymmetric, gate mod- 
ulated and interchangeable, and followed occasionally by 
extended quasiohmic regimes. While our method is par- 
ticularly suited to systems with large charging and small 
coupling, the opposite regime is usually handled pertur- 
batively by SCF-NEGF. Developing the transport for- 
malism for the intermediate coupling regime could be 
nontrivial involving novel physics due to the inter- 
play between charging (localization) and hybridization 
(derealization), and may be crucial to understanding a 
variety of other molecular switching and sensing-based 
phenomena already being explored experimentally. 
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